Finite-size effects on open chaotic advection 
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We study the effects of finite-sizeness on small, neutrally buoyant, spherical particles advected 
by open chaotic flows. We show that, when projected onto configuration space, the advected finite- 
size particles disperse about the unstable manifold of the chaotic saddle that governs the passive 
advection. Using a discrete-time system for the dynamics, we obtain an expression predicting the 
dispersion of the finite-size particles in terms of their Stokes parameter at the onset of the finite- 
sizeness induced dispersion. We test our theory in a system derived from a flow and flnd remarkable 
agreement between our expression and the numerically measured dispersion. 

PACS numbers: 



I. INTRODUCTION 

A proper understanding of chaotic advection [1] in in- 
compressible, open flows is substantially relevant, as ap- 
plications range from laboratory experiments [2] to envi- 
ronmental flows of occanographic [3] and atmospherical 
importance [4]. Usually, open flows displaying chaotic 
advection are asymptotically regular, but the dynamics 
of the advected particles possesses a chaotic saddle [5] 
which governs the transient behavior of the orbits. This 
saddle is an invariant set consisting of infinitely many 
unstable orbits organized in a fractal structure. Its sta- 
ble and unstable manifolds also display fractality and can 
be observed directly in physical space. In particular, the 
unstable manifold is traced out by an ensemble of fluid 
particles initially placed in the region of the saddle. 

The fundamental aspects of chaotic advection in open 
flows are now relatively well known (see [6] and references 
therein) , but almost exclusively in the case where the par- 
ticles are considered to move as fluid particles (passive 
tracers), without inertia. This situation is termed pas- 
sive advection. In many important flows, however, the 
finite-sizeness of the particles has to be considered [7]. 
The resulting dynamics is strikingly different from pas- 
sive advection. In particular, the phase space for finite- 
size particles has twice the number of dimensions of the 
phase space of passive advection. The latter is simply the 
configuration space (physical space), and corresponds to 
an invariant subspace of the former. The higher dimen- 
sionality results from the new degrees of freedom corre- 
sponding to the components of the velocity of the finite- 
size particles. This is a consequence of the fact that the 
finite-size particles are not constrained to have the veloc- 
ity of the advected fluid. By contrast, the velocities of 
passive tracers do coincide with the velocity of the flow. 

One consequence of this difference is that the motion 
of finite-size particles can diverge from the motion of the 
corresponding passive tracers in certain regions of the 
flow, even when the particles have the same density as 
the fluid [8]. We stress that this divergence arises solely 
from the finite-size character of the particles, which we 
consider here to be neutrally buoyant [9], i.e., they have 
the same density as the fluid. The aim of this paper is to 



investigate the consequences of this divergence to open 
chaotic advection. We argue and show that, when pro- 
jected onto the configuration space, the advected finite- 
size particles disperse about the chaotic saddle (and its 
unstable manifold) corresponding to passive advection. 
This is a general phenomenon that had not been reported 
so far. We develop a theory to analyze it. We obtain a 
quantitative expression for the dispersion of the finite-size 
particles about the chaotic saddle and its unstable man- 
ifold in configuration space as a function of the Stokes 
parameter. We test and validate our theory using the 
blinking vortex-source system and find a strikingly good 
agreement with the directly measured dispersion. A ma- 
jor consequence of this dispersion is that finite-size effects 
can destroy the fractal structure of open chaotic advec- 
tion in configuration space, with severe consequences to 
active flows [10, II]. 



II. CONTINUOUS DESCRIPTION OF THE 
DYNAMICS 

In dimensionless form, the equation of motion for a 
small rigid spherical particle advected by an incompress- 
ible fluid of same density with a given velocity field u(r, t) 
is [12] 



dv _ Du 1 1 /dv Du 

di~D^~St^'^""''"2Vrfi "Dt 



(1) 



where v is the velocity of the particle and St = 
is its Stokes number. Here a is the radius of the parti- 
cle, whereas U and L are the characteristic velocity and 
lenght of the flow, respectively. The kinematic viscosity 
of the fluid is given by f. Physically, the Stokes num- 
ber is a measure of the finite-size effects. In the limit of 
vanishing Stokes number, we recover passive advection, 
with V = u. Equation (1) is Newton's law with the terms 
on the right hand side corresponding, respectively, to the 
force exerted by the undisturbed flow. Stokes drag and 
the added-mass effect. In this approximation, the Faxen 
corrections and the Basset-Boussinesq history force term 
are neglected [12]. 
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If we write the derivative along the trajectory of the 
fluid element, ^ = ^ + (u • V)u, in terms of the 
derivative along the trajectory of the particle, ^ = 
1^ + (v • V)u, we cast Eq. (1) into the form [8] 

where A = V — u and Ju is the Jacobian of u. For con- 
venience, we treat a two-dimensional fluid flow, so the 
dynamics of finite-size particles occurs in a 4-dimensional 
phase space. The configuration space, corresponding to 
A = 0, is a 2-dimensional invariant subspace where pas- 
sive advection takes place. 

III. THE OPEN FLOW MODEL 

In order to illustrate the finite-size effects in open 
chaotic advection, we choose the blinking vortex-source 
system for the flow [6, 13]. This system is periodic 
and consists of two alternately open point sources in 
a plane. It models the alternate injection of rotating 
fluid in a large shallow basin. Apart from the two 
point sources, the dynamics is Hamiltonian, with the 
streamfunction given by = — {K \nr' + Qcj)') Q{t) — 
{K\nr" + Q(I}")Q{-t), where r = 0.5T - i mod T. 
Here, r' and (/>' are polar coordinates centered at (—1, 0) 
whereas r" and 4>" are polar coordinates centered at 
(1,0). The two parameters Q and K are, respectively, 
the strenghts of the source and of the vortex. The pe- 
riod of the flow is T, whereas Q stands for the Hcavi- 
sidc step function. The sources are located at positions 
(±1,0). For each half period, the system remains sta- 
tionary with only one of the sources open. This allows 
one to analytically integrate the equations of motion, 
Ux = d^/dy and Uy = —d'i'/dx, during each half pe- 
riod, and thus to write explicitly a stroboscopic map for 
this system, recording the positions of the particles after 
integer multiples of the period T. In complex represen- 
tation, z ^ X + iy = r exp(i(/)), the stroboscopic map 
is 

/ X 1/2-14/2 

^»+i = (4-i)(i-^r^j +1^ (3) 

where z'^ — {zn + 1) {1 — ^ ^^^2 )^/^~'^/^ — 1, and where 
the two parameters governing the dynamics are rj = QT 
and ^ = K/Q. We fix 77 = —0.5 and ^ = 10, for which 
the dynamics is chaotic. Figure 1(a) shows the unstable 
manifold of the chaotic saddle for passive tracers. 
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FIG. 1: (a) Unstable manifold of the chaotic saddle, (b) The 
projection onto configuration space of the seventh iteration 
of Eq. (5) for an ensemble of particles with e ' = 0.2. Ini- 
tially, 10^ particles are uniformly distributed in the hypercube 
[-0.5, 0.5] X [-1, 0] X [-10-^ 10"^] X [-10-^ 10"^] C R*. We 
plot only the region [—2, 2] x [—2, 1], where about 9 x 10'' are 
mapped to. 

by the two-dimensional area preserving map Xn+i 
f(x„), in our case the one given by Eq. (3) with x„ 
written in complex representation Analogously, a 

discrete version of Eq. (2) for the dynamics of finite-size 
particles is given by [14] 

x„+2 - f(x„+i) = e~'^Jf (x„) (x„+i - f(x„)) . (4) 

Equation (4) contains all the essential features of the flow 
given by Eq. (2). In particular, both have constant rate 
of phase space contraction. The term — (Ju + 'J^l) 
Eq. (2) is substituted by the expected term e~'''Jf(x„) 
in the map. The number e^'^ plays the role of the Stokes 
parameter in the discrete dynamics. Equation (4) has 
been successfully used in the study of the dynamics of 
neutrally buoyant finite-size particles under chaotic ad- 
vection [15] and it was generalized in order to describe 
also the dynamics of particles whose density differs from 
that of the fluid [16]. It is convenient to cast Eq. (4) into 

f X„+l = f(x„) + Vif„, 
I W„+i = C ''Jf(x„)v^f„. 

If Ax and A^^ are the eigenvalues of Jf(x), the vec- 
tor w„ = x,i-|_i — f(x„) is amplified in regions where 
7 < I In I Ax II . This corresponds to regions of the flow 
characterized by larger strain. In open chaotic flows with 
asymptotic regularity, such regions are precisely where 
the chaotic saddle is located. Consequently, there the 
finite-size particles may detach from their corresponding 
fluid elements and it is where the finite-size effects are 
expected to occur. 



IV. DISCRETE DESCRIPTION OF THE 
DYNAMICS 

Equation (3) is the discrete version of the flow whose 
streamfunction is Thus, passive advection is described 



V. DISPERSION OF THE ENSEMBLE OF 
FINITE-SIZE PARTICLES 

Analogously to the continuous case, the dynamics de- 
scribed by Eq. (5) takes place in a 4-dimensional phase 



space. The configuration space corresponds to the 2- 
dimensional invariant subspace w = 0, where passive 
advection occurs. The hmit of vanishing particle size 
corresponds to St ^ in the flow and to ^ in 
the map. The solutions in this case are, respectively, 
V = u and x,i+i = f(x,i). In this limit, we recover the 
motion of passive advection and an ensemble of particles 
initially located in the region of the chaotic saddle traces 
asymptotically its unstable manifold in the configuration 
space. As the size of the particles grows from zero, how- 
ever, the particles trace out the unstable manifold in the 
4-dimensional phase space. The key point is that, when 
projected onto the configuration space, they do not trace 
the unstable manifold, W^{T,), corresponding to passive 
advection, but they instead "disperse" about it. Figure 
1(b) shows the projection onto the 2-dimensional config- 
uration space of the seventh iteration of an ensemble of 
finite-size particles, initially placed in the region of the 
chaotic saddle in the 4-dimensional phase space. We thus 
have a clear qualitative picture of the finite-size effects, 
in which the smearing of the fractal structure in config- 
uration space is evident. As the size of the particles gets 
larger, their "dispersion" about increases. 



VI. THEORY FOR THE DISPERSION 

To obtain a quantitative description of the phe- 
nomenon, we define the dispersion D oi a set Si about 
a set 5*2 as the average of the distances d{x, S2) between 
the points x € Si and the set 5*2. So we have 



(a) 



(6) 



where d{x,S2) = iTnii{d{x , y) , y S 6*2}. Our goal is the 
derivation of an expression predicting the behaviour of 
the dispersion, about VF"(S]), of the position vectors in 
configuration space of an ensemble of finite-size particles 
whose initial conditions are (xg, wg), as a function of e""*". 
Differently from the case illustrated in Fig. 1(b), the ini- 
tial positions xq are chosen in the chaotic saddle, S, as 
our aim is to make a theory on the invariant set. Fig- 
ure 2 shows the chaotic saddle. It is carefully obtained to 
warrant the correct natural measure [17]. The initial con- 
ditions Wq are chosen so that ||wo|| ^ 1, corresponding 
to V ~ u in the flow, a condition assumed in the deriva- 
tion of Eq. (1). The directions of the vectors wg are 
randomly chosen with uniform angular distribution. We 
calculate the dispersion for the second iterate of the map, 
as this is the lowest iterate where the flnite-size effects are 
present. The first iterate of the map, xi, is ||wo|| -distant 
from f(xo), independently of e"''', as we can see readily 
from Eq. (5). Figure 3 illustrates the situation that we 
aim to quantify. 

In order to obtain the dispersion D of the ensemble 
of finite-size particles about M^"(S), wc follow two steps. 
First, wc derive an expression for the expected value, over 
both the natural measure of the chaotic saddle and the 
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FIG. 2: (a) The chaotic saddle (black) and its unstable man- 
ifold (grey), (b) Magnification of (a). 



distribution of the initial vectors Wq, of the distance be- 
tween X2 and the unstable subspace arising from f^(xo). 
This expected value depends on the Stokes parame- 
ter e"'''. The derivation of involves only dynamical 
arguments. Second, we derive an expression for D as a 
function of S-y. This expression is necessary because the 
value of D is smaller than 5^. Vl^"(E) is a fractal set and 
usually contains curves whose expected distance to the 
vector positions of the finite-size particles is smaller than 
Sj. The situation is sketched in Fig. 4. Accordingly, 
the correct value of the dispersion D of the finite-size 
particles about M^"(E) is given by the expression 



(7) 



which defines the function F. This function depends only 
on the geometry of the set VF"(E) in the region close to 
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FIG. 3: (a) The projection onto configuration space of the 
second iteration of Eq. (5) for an ensemble of particles with 
e""*" — 0.2. Initially, 6 x 10* particles (black points) are placed 
in the chaotic saddle and the initial vectors wo are uniformly 
chosen in the range < ||wo|| < 10"'', and with uniform 
angular distribution. The unstable manifold is also 

shown (grey points), (b) Magnification of (a). 

the chaotic saddle. While an analytical expression for F 
is a rather nontrivial task, this problem can be overcomed 
by means of a computer-aided approach. 

Let us now derive the expression for S-y. Projected 
onto configuration space, the second iterate of the map, 
to first order in wq, is given by 

X2 = f (xt) + Jf (xt)wo + e-^Jf (xo)wo, (8) 

where x^ ~ f(xo). The term Jf(x^)wo corresponds to a 
point of the ellipse centered at the origin and having axes 
of lenghts equal to -^/AiiHwoH , where Xu, i = 1,2, are the 



eigenvalues of Jf (x*)jj'(x^). In an analogous fashion, the 
term e~'''Jf (xo)wo corresponds to a point of the ellipse 
centered at the origin and having axes whose lenghts are 
e~'*'v'Aoi'||wo|l, where Aoi, i = 1,2, are the eigenvalues 
of Jf (xo)J|'(xo). Thus we see that the right-hand side 
of Eq. (8) corresponds to a sum of two position vectors 
located at ellipses centered at f(x2). 

Now, because of the large strain in the region of the 
chaotic saddle, the average over its natural measure of 
the lenght of the major axis of the ellipse Jf (x)iV, where 
N is the unit disk in R^, is expected to be much larger 
than unit. This fact implies that the major axis of the el- 
lipse ^f{xl)N is, in general, approximately colinear with 
the unstable subspace i?"(f(x|)) = i?"(f^(xo)) arising 
from f^(xo). To understand that, consider the unit disk 
N centered at xj . Let 02 and 62 be the unit vectors in the 
directions, respectivelly, of the major and minor axes of 
the ellipse Jf (x*)iV. Let ai and 61 be the unit vectors in 
the directions of the pre- images of 0,2 and 62, respectively. 
Now let Ui be the unit vector in the direction of the un- 
stable space of xj'. We can write ui = Cidi-|-C2&i, with ci 
and C2 uniquely determined, as ai and bi form a basis in 
the plane. We then have ^^{xl)ul = ci-\/Aa2 + (c2/VA)62, 
where A = max{Aii,Ai2} (notice that An = I/A12, as 
f is an area preserving map). The vector Jf(xJ)Mi is, 
of course, in the same direction of the unstable sub- 
space arising from the point f^(xo). Taking the scalar 
product Jf(x|f)Mi • 0,2, we obtain the angle 9 between 
02 and the unstable subspace arising from f^(xo) as 

^ /2 

9 = arccos [(1 -I- C2/(c^A^)) ]. Because C2 is not ex- 
pected to be much larger than ci and because of the 
large strain in the region of the chaotic saddle, we typi- 
cally have |c2/(ciA)| ^ 1. Thus, to second order in 1/A, 
we have 9 = |c2/(ciA)|, and we see that the major axis 
of the ellipse Jf (x];)iV is usually almost colinear with the 
unstable subspace arising from f^(xo). 

The directions of the axes of the ellipse c~'' Jf{xo)N , on 
the other hand, do not depend on the direction of the un- 
stable space at f2(xo). Fi gure 5 illustrates the situation. 
Assuming that the distribution of the angles between the 
major axes of the ellipses Jf (x*)iV and e^'^'Jf (xo)Af is uni- 
form over the natural measure of the chaotic saddle (our 
theory remains a good approximation as long as this dis- 
tribution is not too concentrated on a small range of an- 
gles) and choosing our reference frame centered at f^(xo) 
and such that the x-axis is in the direction of the unstable 
subspace arising from f^(xo), the expected value for the 
distance between X2 and the unstable subspace arising 
from f^(xo) is given by 

1 />27r p27T fliz 

where r = {x,y) = (A + c^'^Rq • B • R/3).ro, A = 
/1/(A-V2) \ /(AV2) \ 

1^ (A-i/2) - \ l/(Ai/2) j' 

and R/3 are the rotation matrices for the angles a and (3, 
and ro = (||wo||)(cos sinw). The averages {\~'^/'^) and 
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FIG. 4: Sketch of the position X2 in configuration space, and 
its distances di to W^"(E) and d2 to the unstable subspace 
arising from f^(xo). The dispersion D is the average of di, 
whereas is the average of d2- 
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FIG. 5; Sketch of the situation leading to Eq. (9). 



(A^/^) are taken over the natural measure of the chaotic 
saddle. The right-hand side of Eq. (9) represents an av- 
erage over the directions of the position vectors Jf (x*)wo 
(integral over uj) and c~^Jf (xo)wo (integral over and 
over the direction of the major axis of the ellipse ^[{xo)N 
(integral over a) (see Fig. 5). The matrix A accounts for 
the effects of the term Jf (x*)wo. Its contribution to the 
distance from the point X2 to the unstable subspace aris- 
ing from f^(xo) depends on the minor axis of the ellipse 
Jf (xJ)A^. For this reason the non-zero elements of A in- 
volve (A~^/^). The matrix e~''RQ • B • describes the 
effects of the term e~'''Jf (xo)wo. The contribution of this 
term to 6^ depends essentially on the major axis of the 
ellipse Jf (xo)iV. This is why the non-zero elements of B 
involve (Ai/2) (notice that (A^/^) ^ 1/(^-1/2)). 
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FIG. 6: (color online). Theoretical (red x symbols) and 
numerical (black circles) values of the dispersion D of an en- 
semble of 6 X 10* finite-size particles about 1.8 x 10^ points 
representing TV(E). For the theoretical values, we use Eq. 
(7) (along with Eq. (9)). The values (A"^/^) = 0.12 and 
(A^/^) = 21.06 were computed numerically. The inset shows 
the dispersion Z) of 6 x 10'* randomly generated points (at a 
distance 5tt/2 from the points of the chaotic saddle) about the 
same 1.8 x 10^ points representing W^"(S). Interestingly, the 
overall scaling is linear in the range 3 x 10^^ < 5 < 9 x 10^*. 
By regression, we obtain F = 7.128 x 10"^ -f 2.266 x 10~^ x 5. 
The fact that F{0) 7^ is, of course, due to our approximation 
of W^{T,) by a finite number of points. 



points in the chaotic saddle. The factor 7r/2 is necessary 
because, if Stt/2 is the isotropically distributed distance 
to the point at the chaotic saddle, then S is the expected 
value for the distance to the unstable subspace arising 
from that point. 

In order to validate our theory, we also measured nu- 
mcricaly the dispersion, directly using the definition, Eq. 
(6). The results are also presented in Fig. 6, for com- 
parison, and they are very close to the ones predicted by 
our theory. 



VIII. CONCLUSIONS 



VII. RESULTS 

We have computed D using our theory, Eq. (7) (along 
with Eq. (9)). The initial conditions Wq were uniformly 
chosen in the range < ||wo|| < 10"'^, so that (||wo||) = 
5 X 10~^, and with uniform angular distribution. The 
result is shown in Fig. 6. The inset of Fig. 6 regards the 
estimation of the function F. A good approximation for 
F is obtained by assuming that F(S) is the dispersion, 
about M^"(E), of a set of points which are a distance Stt/2 
apart (in homogeneously distributed directions) from the 



In summary, we have derived a theory accounting for 
the quantitative behaviour of the dispersion of neutrally 
buoyant finite-size particles about the unstable manifold 
of the chaotic saddle as a function of the Stokes param- 
eter. Our theory refers to the onset of the finitc-sizeness 
induced dispersion in the discrete description of the dy- 
namics (second iterate of the map given by Eq. (5)). 
The theory involves a dynamical part, where the com- 
putation of the averages (A~^^^) and (A^''^) is required, 
and a geometrical part, regarding the function F. This 
geometrical part will be the same for any theory of the 
dispersion, either for higher iterates of the map, Eq. (5), 
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or for snapshots in the continuous description. Eq. (2). 
It is a highly difficult problem to be solved analitically, 
but it can be overcomed by means of a computer-aided 
approach. The dynamical part of the theory will be dif- 
ferent for higher iterates of the map or for snapshots in 
the continuous description, but it will be based on the 
same type of reasoning. 

An interesting open problem is to determine the re- 
strictions that the dispersion predicted by our theory im- 
poses on the enhancement of activity (e.g. biological) [10] 
of finite-size particles. 
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